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In this work we study the generalization of the problem, considered in [Phys. Rev. E 91 , 013002 
(2015)], to the case of finite correlation time of the environment (velocity) field. The model describes 
a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow. Inertial- 
range asymptotic behavior is studied by means of the field theoretic renormalization group and 
the operator product expansion. The advecting velocity field is Gaussian, with finite correlation 
time and preassigned pair correlation function. Due to the presence of distinguished direction n, 
all the multiloop diagrams in this model are vanish, so that the results obtained are exact. The 
inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite 
correlation time) that correspond to the two nontrivial fixed points of the RG equations. Their 
stability depends on the relation between the exponents in the energy spectrum £ oc k and the 
dispersion law u oc k 2 ff v . In contrast to the well known isotropic Kraichnan’s model, where various 
correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, 
here the corrections to ordinary scaling are polynomials of logarithms of the integral turbulence 
scale L. 
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I. INTRODUCTION 

Over decades much attention has been paid to the 
problem of intermittency and anomalous scaling in fully 
developed turbulence. Both the natural experiments and 
numerical simulations suggest that the violation of the 
classical Kolmogorov-Obukhov theory [l| is even more 
strongly pronounced for a advected Held than for the ve¬ 
locity field itself; see, e.g., i! and references therein. At 
the same time, the problem of passive advection appears 


to be easier tractable theoretically. Although the theoret¬ 
ical description of the fluid turbulence on the basis of the 
stochastic Navier-Stokes (NS) equations remains essen¬ 
tially an open problem, considerable progress has been 
achieved in understanding passive advection by random 
“synthetic” velocity fields. The most remarkable progress 
has been achieved for the so-called Kraichnan’s rapid- 
change model [H, in which the velocity field is modeled 
by Gaussian ensemble, not correlated in time, with zero 
mean and pair correlation function of the form 


(vi(x)vj(x')) = 5{t - t') D 0 f 

Jk> 


Here P^(k) = Sij — kikj/k 2 is the transverse projector, 
k = |k|, Dq > 0 is an amplitude factor, d is the dimen¬ 
sionality of the x space and 0 < £ < 2 is a parameter 
with the real (“Kolmogorov”) value £ = 4/3. For the 
first time, the anomalous exponents have been calculated 
on the basis of a microscopic model and within regular 
expansions in formal small parameters (i3f]. 

A passively advected field may be chosen both scalar 
and vector, the latter case corresponds to the magnetohy¬ 
drodynamic (MHD) turbulence. From the experimental 
point of view it is a special problem, closely related to 
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the processes taken place in solar corona, e.g., with so¬ 
lar wind; for detailed discussion see GHzi and references 
therein. 

In solar flares, highly energetic and anisotropic large- 
scale motions coexist with small-scale coherent struc¬ 
tures, finally responsible for the dissipation. A simplified 
description of the situation was proposed in [6j: the large- 
scale field Bf = riiB 0 dominates the dynamics in the 
distinguished direction n, while the activity in the per¬ 
pendicular plane is described as nearly two-dimensional. 

The observations and simulations show that the scal¬ 
ing behavior in the solar wind is closer to the anomalous 
scaling of the three-dimensional fully developed hydro- 
dynamic turbulence, rather than to simple Iroshnikov- 
Kraichnan scaling suggested by the two-dimensional pic- 
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ture with the inverse energy cascade Q. Thus, further 
analysis of more realistic three-dimensional models is wel¬ 
come. 

One of the possibilities to make original Kraichnan’s 
model tnj anisotropic is to replace the ordinary trans¬ 
verse projector with the tensor quantity Tij( k), which 
contains a fixed unit vector n: 

Tij( k) = a(ip)Pij{ k) + b(ip)n s niPi s (k)Pji(k), (1.2) 

where a(if) and b(if) are some functions of jb, the angle 
between the vectors n and k; see, e.g., |8l4l0l|. This for¬ 
mulation of the problem corresponds to the small-scale 
anisotropy and contains an isotropic model as a special 
case, if = 1 and b(ip) = 0. 

Another possibility is the “strongly anisotropic” model 
that does not contain an isotropic one as a special case 
and is obtained by introducing the velocity field v having 
preferred direction n: 

»(i,x)=nx«(l, xj_). (1-3) 

In this paper, we consider a more realistic model with 
finite (and not small) correlation time. For this purpose 
the correlation function m has to be modified, and in¬ 
stead of a constant, which is Fourier transform of 6(t—t '), 
in the frequency space it becomes a function of w. In com¬ 
mon cases this modification disrupts the Galilean invari¬ 
ance EH and is interesting only as a model, but in the 
presence of the anisotropy Galilean invariance survives 
and the model is invariant under some special Galilean 
transformations (more precisely see below). 

The energy spectrum of the velocity in the inertial 
range has the form £ oc while the correlation time 

at the momentum k scales as k~ 2+v . Such ensemble 
was employed in some models, studied in [l2], [k|- It 
was shown that, depending on the values of the expo¬ 
nents £ and r), the model reveals various types of inertial- 
range scaling regimes with nontrivial anomalous expo¬ 
nents, which were explicitly derived to the first EH and 
second [l3[ orders of the double expansion in £ and rj. 

It is necessity to stress, that the Kraichnan’s 
model (O) and its generalizations correspond to pas¬ 
sive field approximation: if we neglect the influence of 
advected field 6 to the dynamics of the environment (ve¬ 
locity) field v , the latter can be modeled by statistical 
ensembles with prescribed properties. This approxima¬ 
tion is valid when the gradients of the magnetic fields are 
not too large. 

A most powerful method to study the anomalous scal¬ 
ing in various statistical models of turbulent advection 
provided by the field theoretic renormalization group 
(RG) and operator product expansion (OPE); see the 
monographs El Il5| and references therein. In the 
RG+OPE scenario [3, anomalous scaling emerges as 
a consequence of the existence in the model of compos¬ 
ite fields (“composite operators” in the quantum-field 
terminology) with negative scaling dimensions; see EH 
for a review and the references. In a number of pa¬ 
pers the RG+OPE approach was applied to the case of 


passive vector (magnetic) fields in Kraichnan’s ensemble, 
and to its generalizations (large-scale anisotropy, helicity, 
compressibility, finite correlation time, non-Gaussianity, 
more general form of the nonlinearity); see fl8l - [22| and 
references therein. Explicit analytical expressions were 
derived for the anomalous exponents to the first ITbI] and 
the second [lf| [2(1 orders in £. For the pair correlation 
function of the magnetic field, exact results were obtained 
within the zero-mode approach |23j |. 

In this paper, we apply the field theoretic renormal¬ 
ization group and operator product expansion to the 
inertial-range behavior of strongly anisotropic MHD tur¬ 
bulence within the framework of a simplified model, 
which corresponds to the problem of a passive vector 
field advected by the Gaussian ensemble with prescribed 
statistics. The velocity field v is chosen to be oriented 
along a fixed direction n (“orientation of a large-scale 
flare” in the context of the solar corona dynamics) and 
depends only on the coordinates in the subspace orthogo¬ 
nal to n. In the momentum space, its correlation function 
is some function of k± and frequency w, where k± = |kx| 
and kj_ is the component of the momentum (wave num¬ 
ber) k perpendicular to n. This model can be viewed as 
a d-dimensional generalization of the strongly anisotropic 
velocity ensemble introduced in 24] in connection with 
the turbulent diffusion problem and further studied and 
generalized in a number of papers [251 - 129) . 

The advecting equation for the passive field 6 involves 
a general relative coefficient A , which unifies different 
physical situations: the kinematic MHD model, the lin¬ 
earized NS equation and the passive admixture with com¬ 
plex internal structure of the particles. 


In [2f| the problem of anomalous scaling in the higher- 
order correlation functions of a scalar field, advected by 
such a velocity ensemble, was studied by the RG+OPE 
techniques. It was shown that there exists some set of 
fixed points, which governs infrared (IR) behavior of the 
system. Another conclusion of that work is that in sharp 
contrast to the isotropic Kraichnan’s model and its nu¬ 
merous descendants, due to the mixing of families of rel¬ 
evant composite operators the correlation functions show 
no anomalous scaling and have finite limits when the in¬ 
tegral turbulence scale tends to infinity. 

Further modification of that problem, namely advec¬ 
tion of the vector field by decorrelated in time velocity 
field, was studied in [3Cj. In contrast to [29[ , the inertial- 
range behavior of vector fields appears to be even more 
exotic: instead of power-like anomalies, there are log¬ 
arithmic corrections to ordinary scaling, determined by 
naive (canonical) dimensions. 

The main result of the present paper is that the 
inertial-range behavior of vector fields advected by veloc¬ 
ity ensemble with finite correlation time combines both 
the above features: as in the scalar case, there is a set 
of fixed points, governing the IR behavior; as in the 
zero-time correlation model, the inertial-range behavior 
of vector fields has logarithmic corrections to ordinary 
scaling. The key point is that the matrices of scaling 
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dimensions (“critical dimensions” in the terminology of 
the theory of critical state) of the relevant families of 
composite operators appear nilpotent and cannot be di¬ 
agonalized. They can only be brought to Jordan form; 
hence the logarithms. 

Another interesting property, inherited from the zero¬ 
time correlation model, is that all multiloop diagrams 
are equal to zero and therefore the set of fixed points 
and the existence of logarithmic corrections are proven 
exactly. Moreover, in contrast to previous one, this model 
has two types of such nontrivial diagrams, with different 
causes to be equal to zero. The physical meaning of this 
feature is not yet clarified, but it is clear that it is closely 
connected with the presence of the anisotropy vector n. 

The paper is organized as follows. 

In Sec. eh we give a detailed description of the model. 
In Sec. lIIIl we present the field theoretic formulation of the 
model and the corresponding diagrammatic techniques. 
In Sec. lIVI we establish renormalizability of the model and 
derive explicit exact expressions for the renormalization 
constants and RG functions (anomalous dimensions and 
/3-functions). Due to the presence of the anisotropy, the 
linear response function, the only Green function in the 
model that contains superficial ultraviolet (UV) diver¬ 
gences, is given exactly by the one-loop approximation. 

It is shown that the IR behavior of the model is con¬ 
fined with only two limiting cases: the rapid-change type 
behavior and the “frozen” (time-independent) behavior. 
In contrast to the isotropic case, where the physical 
(Kolmogorov) point £ = 8/3, rj = 4/3 lies exactly on 
the crossover line between the rapid-change and frozen 
regimes [l^, IU, S3 > now this point lies deep inside the 
domain of stability of the nontrivial rapid-change behav¬ 
ior; there is no crossover line going through this point. 
This result is in agreement with the exact analysis of the 
d = (1 + l)-dimensional case [13 and in disagreement 
with [24|, [25| . 

The corresponding differential equations of IR scaling 
are derived, with the exactly known critical dimensions. 

In Sec. El we discuss the renormalization of composite 
operators and present explicit expressions for the matri¬ 
ces of anomalous dimensions and critical dimensions. It 
is shown that these matrices are given exactly by the 
one-loop approximation. The matrices of anomalous di¬ 
mensions appear to be nilpotent. As a result, the IR be¬ 
havior of the pair correlation functions of the composite 
operators is given by canonical powers, corrected by poly¬ 
nomials of logarithms. To obtain inertial-range behavior 
we have to combine this result with the corresponding 
OPE’s. Finally, asymptotic behavior of the pair corre¬ 
lation functions involves two types of large logarithms, 
where the separation enters with the typical UV and IR 
scales (dissipation scale and integral scale). 

Sec. EH is reserved for conclusions. 


II. DESCRIPTION OF THE MODEL 

If the field v is chosen in the strongly anisotropic form 
OD, the turbulent advection of a passive vector field 
Q[x ) = 0(t, x) is described by the stochastic equation [301, 

ia 

d t 0i + d k ( v k 0i - A 0 ViO k ) + dV = v 0 + fod^)6 z + f u 

( 2 . 1 ) 

where 9i(x ) is a vector field, x = {t,x}, d t = d/dt, 
di = d/dxi, n is a unit vector that determines the dis¬ 
tinguished direction, x_l and d± are the components of 
the vectors x and d perpendicular to n, <9|| = d n, v 0 
is the molecular diffusivity coefficient, d 2 is the Laplace 
operator, v(x) = {ui(:c)} is the velocity field, ./) = fi{x) 
is an artificial Gaussian scalar noise with zero mean and 
correlation function 

x) f k (t', x')) = S(t — t') C ik (v/L). (2.2) 

Here r = x — x', r = |r|, the parameter L = M~ x is the 
integral (external) turbulence scale related to the stirring, 
and Ci k is a dimensionless function finite for r/L —> 0 and 
rapidly decaying for r/L —>• oo. 

Both v and 6 are divergence-free (“solenoidal”) vector 
fields: 

diVi = 0, diOi = 0. (2.3) 

Following j33j, we included into the stochastic 
advection-diffusion equation ( 12 . 11 ) additional arbitrary di¬ 
mensionless parameter Mo, which unifies different physi¬ 
cal situations: the case Mo = 1 corresponds to the kine¬ 
matic MHD equation, describing, for example, the evo¬ 
lution of the fluctuating part 6 = 6{x) of the magnetic 
field in the presence of a mean component 9 °, which is 
supposed to be varying on a very large scale; the case 
Mo = — 1 corresponds to the linearization of the NS equa¬ 
tion around the rapid-change background velocity field; 
in the case Mo = 0 equation ( 12 . 11 ) loses the stretching 
term d k {yi6 k ) and the model acquires additional symme¬ 
try under translations 9^9 + const. This case has to 
be studied separately, see j34|. 

The pressure term can be expressed as the solution of 
the Poisson equation 

d 2 V = (Mo - 1) div k d k 9i (2.4) 

and is needed to reconcile dynamics of the field 6i with 
transversality condition & 

For renormalizability reasons it is necessary to intro¬ 
duce additional dimensionless constant /o, which breaks 
the Od symmetry of the Laplace operator to Od-i <8> Z 2 : 
d 2 —> d\ + /o<9j| (Z 2 is the reflection symmetry x\\ —» 
—£|i). Interpretation of the splitting of the Laplacian 
term can be twofold; cf. [29[. On one hand, stochastic 
models of the type (12.11) must include all the IR relevant 
terms allowed by the symmetry, therefore it is natural to 
include the general value /o ^ 1 to the model from the 
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very beginning. On the other hand, the extension of the 
model to the case / 0 7 ^ 1 can be viewed as a purely techni¬ 
cal trick which is only needed to ensure the multiplicative 
renormalizability and to derive the RG equations. 

Instead of the real problem, where the velocity field 
v(x) has to satisfy the NS equation with some additional 
terms that describe the feedback of the advected field 
9(x) on the velocity field, we will consider the kinematic 
problem, where the reaction of the field 9(x) on the ve¬ 
locity field v{x) is neglected. It is assumed that, if the 
gradients of 9(x) are not too large, it does not affect es¬ 
sentially dynamics of the conducting fluid. Thus, the field 
v(x) can be simulated by statistical ensemble with pre¬ 
scribed statistics. It is assumed to be Gaussian, strongly 
anisotropic [see EDI, homogeneous, with zero mean and 
a correlation function Qjj, [u| |29| 


(Vi(t, X) V k (t', x')) = UiHk (v(t, X_j_ ) v(t’, x'jJ), 

(2.5) 

where 

(v(t, xj.) v(t', Xj_)} = f e * k -( x - x ') D v (ui, k). 

(2-6) 

The function D v is chosen in the form 


D v (ix,k) = 27rd(fc||) Do 


'.5 -d-(£+n) 


a 0 vok 


2-r, 

± 


1 2 • 


(2.7) 


Here d is the dimensionality of the x space, k± = |kj_|, 

I 


5($) 


■ 2 V iD v X Vk + -O^DgO'k + 6' k 


-d t e k 


Here all the terms, with the exception of the first one, rep¬ 
resent the De Dominicis-Janssen action for the stochastic 
problem (ED , ED at fixed v , while the first term rep¬ 
resents the Gaussian averaging over v. Furthermore, Dg 
and D v are the correlators ED and (ED respectively; 
the needed integrations over x = (t, x) and summations 
over the vector indices are implied. 

As a rule, synthetic velocity ensembles with finite cor¬ 
relation time suffer from the lack of Galilean invariance, 
which can lead to some physical pathologies; see, e.g., 
the discussion in 0 - Surprisingly enough, the presence 
of the anisotropy can improve the situation. 

Indeed, it is directly checked that in our strongly 
anisotropic case the action functional ED with the cor¬ 
relator ED in its first term appears invariant with re¬ 
spect to the Galilean transformation of a special form: 

6(t, x) —» 0(t, x + ut), 0'{t , x) —>• 6'{t, x + ut ), 

v(t, x) — > v(t,K + ut) — u. (3.2) 


1 /m is another integral turbulence scale, related to the 
stirring, D 0 > 0 is an amplitude factor and symbol fcy 
denotes the scalar product k • n. The function (12.611 in¬ 
volves two independent exponents £ and 77, which in the 
RG approach play the role of two formal expansion pa¬ 
rameters; a new parameter ao is needed for the dimen¬ 
sionality reason. Depending of this parameter, the func¬ 
tion ED demonstrates two interesting limiting cases: if 
ao —>■ 0, D v (lo) oc d(ui), so that from the physs point 
of view this situation corresponds to the independent of 
time (“frozen”) velocity field. The situation ao —t 00 in 
fact means that (ao^o ) 2 w 2 , so that this case corre¬ 

sponds to the rapid-change model. 

The relations 


Do/vlh = go = A« + " (2.8) 

define the coupling constant go, which plays the role of 
the expansion parameter in the ordinary perturbation 
theory, and the characteristic UV momentum scale A. 

III. FIELD THEORETIC FORMULATION OF 
THE MODEL 

A. The action functional and the Galilean 
symmetry 

The stochastic problem ED ED is equivalent to the 
field theoretic model of the extended set of three fields 
<I> = {9,9', v} with the action functional 


- ( Vidi)0 k + Ao(9idi)v k + u 0 (dj_ + fodf\)9 k 


(3.1) 


Here the transformation parameter has the form u = n u 
with the vector n from in, so that the scalar coeffi¬ 
cient in (11.31) changes as v(t, x_l) —> v(t , xj_) — u and the 
arguments xy of all the fields in (13.21) remain intact. 

This fact can be interpreted as follows. Consider the 
generalized stochastic NS equation 


d t Vi + (yidi)vi + dtp = Rvi + fa, (3.3) 


where R is some differential operation acting only on spa¬ 
tial coordinates and p = —d~ 2 (diVi)(diVi) is the pressure. 
If the random force fa is taken to be white in time, the 
equation ED is Galilean covariant because it involves 
the full covariant derivative d t + (■ vidi ). 

However, for the velocity field of the form (11.31) all the 
nonlinear terms in (13.31) vanish due to the independence 
of the scalar coefficient v on ary: v k d k Vi = ntvd\\v = 0, 
and similarly for the pressure. Thus the equation (13.31) 
becomes in fact linear and generates a Gaussian velocity 










5 


field. Its pair correlation function has the form 


(v*Vj) 


fa 2 + R 2 ( k) ’ 


(3.4) 


where D)’ (k) is the pair correlator of the random force 
tfa. It coincides with (12.51) if one choses (in the momen¬ 
tum representation) I?(k) = uoVok\ v and fa = 4>rii with 

{(/></)) = ffo^o/o S(t — t') S(k^)k^ d ~ < ' E+r> \ It remains to 
note that the resulting velocity ensemble has a finite cor¬ 
relation time in contrast to the random force fa in (13.31) . 


Here CW(k) oc Pw{ k) is the Fourier transform of the 
function from EH ; the propagator {0' d 6' d ,) is equal to 
zero. 

In fact, the action functional EH has to be modified 
for the sake of renormalizability. As a consequence, the 
functions m and EH will acquire certain additional 
terms. However, it turns out that those additional terms 
do not contribute to the divergent parts of all the relevant 
diagrams, and thus they can be neglected. These issues 
are discussed in detail in sec. IIVD1 and in the following 
we will use for the propagators the above expressions 
EH) and (10) . 


B. Feynman diagrammatic technique 


C. Canonical dimensions and UV divergences 


The model EH corresponds to a standard Feyn¬ 
man diagrammatic technique with the triple vertex 
d' [—(vidi)6k + Ao{Oidi)vk] and the three bare propaga¬ 
tors. A fragment of arbitrary diagram is represented in 

Fig. (HJ. 



FIG. 1. The triple vertex with three attached propagators. 


The analysis of UV divergences is based on the anal¬ 
ysis of canonical dimensions of the 1-irreducible Green 
functions. In general, dynamic models have two scales: 
canonical dimension of some quantity F (a field or a pa¬ 
rameter in the action functional) is completely character¬ 
ized by two numbers, the frequency dimension dp and the 
momentum dimension dp. They are defined such that 

[F]~[T]- d F[L]- d F, (3.8) 

where L is some reference length scale and T is a time 
scale. 

In the scalar version of strongly anisotropic 
model m eh , however, there are two independent 
length scales, related to the directions perpendicular 
and parallel to the vector n [29|. But the transversality 
conditions 


In the frequency-momentum representation the triple 
vertex corresponds to the expression 


Vcab — k a i“Ao5 a c k^ , ( 3 . 5 ) 


where k s is the momentum of the field 6'\ in the di¬ 
agrams it is represented by the point, in which three 
lines connect with each other. The three propagators 
are determined by the quadratic (free) part of the ac¬ 
tion functional and are represented in the diagrams as 
slashed straight (the slashed end corresponds to the field 
6 '), straight (the end without a slash corresponds to the 
field 9) and wavy (which corresponds to the field v) lines, 
respectively; cf. [3(l |. 

The line ( v a v a ') 0 in the diagrams corresponds to the 
correlation function EH, and the other two propaga¬ 
tors in the frequency-momentum representation have the 
forms 


(0c9' c ,) 0 = 

Pcc'(k) 

(3.6) 

-iuj + vq (k)[ + fok^j 

( 9bOb') 0 = 

C w { k) 

r / \ n 2 ' 

(3.7) 


U) 


( k l + /o*f) 


diOi = 0, = 0 (3.9) 

forbid this option; see [30l |. In particular, this means that, 
in contrast to the scalar case, the constant fo from (12.11) 
in our case is dimensionless. 

The dimensions in (| 3.8 D are found from the obvious 
normalization conditions d\ = —d^ = 1, = — g?£ = 0, 

= —d% = 1, d* = df =0, and from the requirement 
that each term of the action functional EH be dimen¬ 
sionless (with respect to the two independent dimensions 
separately). Based on dp and dp, one can introduce 
the total canonical dimension dp = dp + 2 dp (in the 
free theory, d t oc d\ oc <9y), which plays in the theory 
of renormalization of dynamic models the same role as 
the conventional (momentum) dimension does in static 
problems; see, e.g., [u|. 

The canonical dimensions for the model (13.11) are given 
in Table HI including renormalized parameters, which 
will be introduced a bit later. From Table Q] it follows 
that our model is logarithmic (the coupling constants 
go ~ and op ~ are dimensionless) at 

£ = rj = 0, so that the UV divergences in the Green 
functions manifest themselves as poles in £, 77 and their 
linear combinations. 
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TABLE I. Canonical dimensions of the fields and parameters 


F 

o' 

0 

V 

M, m, /q A 

v,v 0 

A, Ao 

fjo 

ti, Uo 

ao 

go, go 

a, g, g 


1/2 

- 1/2 

1 

0 

1 

0 

0 

0 

0 

0 

0 

~df 

d 

0 

-1 

1 

-2 

0 

0 

0 

V 


0 


d + 1 

-1 

1 

1 

0 

0 

0 

0 

V 

t + v 

0 


The total canonical dimension of an arbitrary 1- 
irreducible Green function Tjv # = (<!>... <1>) i_i r is given 
by the relation 

dr Nit = d+ 2 — 'y ' = d-\- 2 — Ngidgi — Ngdg — N v d v . 

(3.10) 

Here Ng, = { Ng , Ng*, N v } are the numbers of corre¬ 
sponding fields entering the function T, and the sum¬ 
mation over all types of the fields in (13.101) and analogous 
formulae below is always implied. 

Superficial UV divergences, whose removal requires 
counterterms, can be present only in those functions T 
for which the “formal index of divergence” dr Nit is a non¬ 
negative integer. Dimensional analysis should be aug¬ 
mented by the following considerations: 

(1) In any dynamical model of type (13.11) . 1-irreducible 
diagrams with Ng* = 0 necessarily contain closed circuits 
of retarded propagators (13.61) or at least one vanishing 
propagator (9(9^) and therefore vanish. 

(2) For any 1-irreducible Green function Ng> — Ng = 
2 Ad, where Ad > 0 is the total number of the bare prop¬ 
agators (66 )o entering into any of its diagrams. 

(3) Using the transversality condition of the fields 

6i and Vi we can move one derivative from the vertex 
—8' k (vidi)6k + Aq 9' k (8idi)vk onto the field 8(. Therefore, 
in any 1 -irreducible diagram it is always possible to move 
the derivative onto external “tail” 8 ' k , which reduces the 
real index of divergence: d' rN = — Ng>. The field 6' k 

enters the counterterms only in the form of the derivative 
di0' k - 

From Table |H and (13.101) we find that 

dr t ,.=d + 2-{d + \)N e >+N e -N v (3.11) 


Another possibility is (O'6 . .. 9v ... u)i_i r with Ngi = 
1 and arbitrary Ng = N v , for which dr = 1, d' T = 0. 
From the requirement Ng > Ng* it follows that the only 
nonvanishing function of this type is (9' a 9pVy)\.i T . 

IV. RENORMALIZATION OF THE MODEL 

A. Perturbation expansion for the 1-irreducible 
linear response function 

The field theoretic formulation means that statistical 
averages of random quantities in the stochastic prob¬ 
lem ED, ED coincide with functional averages with 
weight exp<S(<E>) with the action (13.11) . 

Let us denote the generating functional of the normal¬ 
ized full Green functions G = ($...$) as G(A ), where 
A(x) = {A(x),A'(x),A v (x)} is the set of “sources,” ar¬ 
bitrary functional arguments of the same nature as the 
corresponding fields. Thus, the generating functional of 
the 1-irreducible Green functions is obtained using the 
Legendre transform: 

r($) = lnG(A)-$A; (4.1) 

see, e.g., [ 3 . 

The Green functions with the auxiliary field O' rep¬ 
resent, in the field theoretic formulation, the response 
functions of the original stochastic problem, in particu¬ 
lar, the simplest (linear) response function is given by the 
relation 

(S0p/6f a ) = (8p & a ). (4.2) 

Let us consider the 1-irreducible linear response func¬ 
tion 


and 


d!r Nt) = {d + 2){l-Ng,) + No-N v . (3.12) 

From these expressions we conclude that, for any d, 
superficial divergences can be present only in the 1 - 
irreducible functions of two types. 

The first example is provided by the infinite family of 
functions (O’ 9... 9) i_j r with Ng* = 1 and arbitrary Ng, 
for which dr = 2, dp = 0. However, all the functions 
with Ng > Ng / vanish (see above) and obviously do not 
require counterterms. Therefore the only nonvanishing 
function from this family is (9' a 9f 3 )i_i r . 


(4.3) 

3>=0 

In accordance with ED generating function for it con¬ 
sists of two parts, 


na/3 


S 6 


2 - (Wl-ir - 5e , a se m) 


r($) = s($) + r($), (4.4) 

where for the functional arguments we have used the 
same symbols = {6,0',v} as for the corresponding 
random fields; 5(4>) is the action functional (13.11) and 
r(4>) is the sum of all the 1 -irreducible diagrams with 
loops. Thus, one obtains 
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= iwP a p{ p) - voP 2 ±P a p(p) - i/ 0 / 0 (pn) 2 P a( g(p) + Y, a p, (4.5) 


where P a p{ p) = d a) g — PaPp/p 2 is transverse projector 
and £ aj g is the “self-energy operator,” diagrammatic rep¬ 
resentation for which is represented in the Fig. [2] Here 



FIG. 2. Diagrammatic representation for £<*£. 


the ellipsis stands for the 2-, 3- and other N-loop dia¬ 
grams. 

The typical feature of all rapid-change models like ED 
with retarded bare propagator of the type (13.611 is that 
all the skeleton multiloop diagrams entering into the self¬ 
energy operator contain closed circuits of such retarded 
propagators and therefore vanish [01 Gil [30|. The depen¬ 
dence of the frequency in function D v [see m] destroys 
this easy construction, and now all the N-loop diagrams 
are expected to give some nontrivial contribution to the 
function T, a p. 

Let us start with the one-loop diagram. It is repre¬ 
sented by the expression 


J a(3 


-*/£. 


dk 


27t S(k\\)k 


5-d-(5+i?) 

X 


(27r)° 


(-K*> + v 0 (p + k)^_ + fo (p + k )jj ) 


U) 2 + 


a 0 v 0 k\ n 


■P a i(p)JijPjp{p), (4.6) 


where the fraction is a product of the propagator func¬ 
tion m and the correlator m , transverse projectors 
P a i( p) and Pgj (p) are present due to the transversal- 
ity conditions (13.91) . and J t] is an index structure of this 
diagram: 


Jij =Viab(p)V dc j(p + k)P bd (p + k)n a n c . (4.7) 

Here and below Vjjfc(p) is the triple vertex (13.51) : the 
Greek letters a, /3 and the Roman letters a-d denote 
the vector indices of the propagators (12.51) and (13.61) with 
the implied summation over repeated indices. Since the 
index of divergence for this diagram dr = 2, we need to 
calculate only the terms, proportional to p 2 . 

The calculation of this diagram is similar to the zero¬ 
time correlation case [3(|, so we will discuss it here only 
briefly. 

The integration over the frequency ui is trivial. In order 

I 


^cn/3 


ffo^o/o 

2q?o 


C d -1 


d — 2 + Aq 


P a p(p) + 


where C d - 1 = S d - i/(27r) d 1 and the vector hk, which is 
orthogonal to p, is defined as 

h k = Pmk{p)n m = n k ~P\\Pk/p 2 - (4.11) 

The integral over k± in expression (14.101) can be sim¬ 
plified in the minimal subtraction (MS) renormalization 


to integrate over the vector k with the function d(/c||) in 
the integrand we need to average the expression (1131) 
over the angles: 


I dk<5(fc||)/(k) 


S d -i / dk ± k d ± ~ 2 (/(k ± )), (4.8) 


where (• • •) is the averaging over the unit sphere in the 
(d — l)-dimensional space, S d - i is its surface area, and 
k± = |kj_ |. To average some function of k± over the 
angles in the orthogonal subspace we use the following 
expression: 


k t k i \ = PiM) 
k\ j (d-1)- 


This gives: 


(4.9) 


Mo — 1) 
d- 1 


2 


(pn) 2 



dk_i 


k 


i-£ 


k 2 , + «o k 


(4.10) 


scheme, which we adopt in what follows. In that scheme, 
all the anomalous dimensions 7 are independent of the 
regularizators like £ and ip and we may chose them arbi¬ 
trary with the only restriction - our diagrams have to re¬ 
main UV finite; see [L§] for detailed discussion. The most 
convenient way is to put 77 = 0, so the expression (14.101) 

























turns into 


v _ ffo^o/o 7 , d — 2 + A t 
2 a 0 (l + a 0 ) d 1 d— 1 


o D ,• A , (A — l) 2 ~ ~ 1 , ^2 f°° 1 

A/3(pH -—n Q n^J (pn) J 




and we obtain the following result: 

ffo^o/o 


Ziap — 


2 ao(l + «o) 


C d _i 


d — 2 + A (A - l) 2 - - 

P a p(p) H-^—n a n/3 


d — 1 


(P ' n ) 2 


7-« 


(4.12) 


(4.13) 


The remaining multiloop diagrams will be discussed a 
bit later, in section HV Cl 


FIG. 3. Diagrammatic representation for A a ^ 7 . 


B. Perturbation expansion for the 1-irreducible 
function (9' a 9 pv-,} ±.i r 

The expansion like (14.51) for the function (9' a 9 pv^) i. n - 
has the form 

(9 a 9p'Vry') I-Iy — Va/3 7 T A a 0^ 

— iS a *yP[3 ^>^0^q/3P7 T f 377 (d-14) 

where V a p^ is the vertex (13.51) and A a p 1 is represented 
in the Fig. [3] 



As in the case of self energy operator in Fig. [ 2 J the 
ellipsis stands for the 2-, 3- and other N-loop diagrams. 

Since our model is Galilean invariant, as discussed in 
Sec. M the terms 9 k d t 9 k and 9' k (vidi)9 k in the action 
functional may be renormalized only with the only renor¬ 
malization constant Z\. The index of divergence for this 
function is dp = 1 , so that the counterterms with d t 
are forbidden. Consequently, counterterm 9 k (vidi)9 k is 
also forbidden. If Ao = 1, the vertex (13.51) is transverse, 
the nonlocal term dV in the stochastic equation (ED is 
absent and the action functional is local in-time. This 
means that the counterterm 9' k (9ndi)v k is forbidden be¬ 
cause the appearance of some constant Z 2 , which this 
term is renormalized by, is equivalent to appearance of 
some multiplier like A 7 ^ 1 , i-e., the appearance of non¬ 
local terms in the action functional. Similar reasoning 
exclude the appearance of such a counterterm if Aq = 0 . 
Thus, we may conclude, that A a ^ 7 is proportional to 
A(A — 1) and vanish for the aforementioned cases. 

The procedure of calculating the one-loop approxima¬ 
tion of A„ p 1 is similar to the one-loop contribution to the 
self-energy operator T, a p, discussed in previous section. 
The analytical expression for the former is 


A a ^3 7 — Dq 


du 

27r 


dk 


( 27r ) d (-iu + v 0 (k + q)^ + / 0 (fc + g)| + u 0 (k - p)^_ + f 0 (k - p)jj ) 


27T<5(fc||)£t L 


5-d-^+p) 


u 2 + 


a 0 v 0 k 


,2-v 


n 2 


Pai (q) Pj'y (p)^t/3, (4.15) 


where p and q are two external momenta, Jipj is the 
index structure of this diagram, transverse projectors 
Pi a { p) and Pjryip) and vector np are present due to the 
transversality conditions (13.91) and definition (11.31) . Since 
the index of divergence for this function is dr = 1 , we 


need to calculate only the term, proportional to the lin¬ 
ear combination of p and q. Also we may put 77 = 0 in 
this diagram and left with the only regularizator £. 

The integral over w is convergent; direct calculation 
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shows that 

J iP j oc A)(l - A 0 )niPpj(n). (4.16) 

This means that 

Ja/3"f = Pai (q) JipjPj'y (p)tt/3 — 0; (4.17) 

i.e., the function (9' a 9 f 3 i> 7 )i_i r does not diverge not only 
for the cases Ao = 0 and Ao = 1, discussed above, but 
also in all the other situations. 

The multiloop diagrams will be discussed in the next 
subsection. 


C. Multiloop diagrams 


Since in any 1-irreducible diagram it is always possible 
to move the derivative onto external “tail” 9 ' k , the real 



FIG. 5. One of two possible fragments of arbitrary multiloop 
diagram for self-energy operator T, a p. 



In order to renormalize our model we have to deal 
with two types of multiloop diagrams - one of types cor¬ 
responds to the function (9' a 9 t 3 )i_i r and is represented 
in Fig. [2] the other one corresponds to the function 
{6' a 9pVry) i_i r and is represented in Fig. [3l Let us start with 
the latter. Any multiloop diagram of this type contains 
a part with the structure, represented in Fig. QJ Since it 


a Y 


b 



c 


a 


FIG. 4. Fragment of arbitrary multiloop diagram, entering 
into expansion of the function (9' a 6pvA\-ir- 

is sufficient to calculate all the diagrams at external mo¬ 
menta equal to zero (the real index of divergence d' r = 0 ), 
the integral, corresponding to the divergent part of the 
diagram, necessarily contains as a factor the following 
expression: 

I 0 oc ( 5 (fc||)( 5 (g||)n a I 4 ac (k)n a V 3 a7 (k + q)P 76 (k), (4.18) 

where V is the vertex (13.51) . and the (5-functions appear 
from velocity correlator m - Since Iq is proportional to 
the sum of k\\ and q\\ with some coefficients, after inte¬ 
gration with the <5-functions all these diagrams vanish. 

Any multilop diagram, entering into the expansion of 
the 1-irreducible linear response function (9' a 0p) i_i r , con¬ 
tains a part with structure, represented in Fig.[5]or a part 
with structure, represented in Fig. [ 6 ] 


FIG. 6. Another possible fragment of arbitrary multiloop di¬ 
agram for self-energy operator Y, a p. 


index of divergence for this diagram d' T = 1. This means, 
that in course of calculation of the structures, represented 
in Fig.[5]and Fig. [ 6 ] we are interested only in terms, linear 
in the external momenta p. 

The analytical expression for the first structure, de¬ 
noted by / 1 , is proportional to 

h oc<5(/c||)<5(<7||)n 

zVxyz (p + k)P yp (p + k - q) (4.19) 
x n q V pqr (p + k - q)P rt (p - q)n n I4„ m (p - q). 

Here p is the external momentum, k and q are inter¬ 
nal integration momenta, V is the vertex (13.51) . P is 
the transverse projector, and the unit vector n and 5- 
functions stem from velocity correlator (12.51) . 

Direct calculation shows, that p is proportional to 
some linear combination of fc|| and , and, as well as in 
the case of Iq, after the integration with the (5-functions 
all diagrams with this structure vanish. 

Another structure, represented in Fig. [ 6 l possess the 
same property - analytical expression for it is similar 
to (14.191) , and, as can be seen from the direct calculation, 
all the diagrams with this structure also appear to be 
equal to zero. 

It should be stressed that, in contrast to rapid-change 
models like m with <5 functions in time, where all these 
multiloop diagrams vanish due to the closed circuits of 
retarded propagators, in our model their vanishing has a 
rather nontrivial origin and results from the presence of 
the anisotropy in it. 


D. Renormalization and RG equations 

Substitution of the explicit expression (14.131) for the 
divergent part of the self-energy operator Y, a p into the 
expression (14.51) for the 1-irreducible linear response func¬ 
tion rf gives 
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-iCk/3 


= {iu-v 0 p 2 ± -i/ 0 fo(p • n) 2 } P a p(p)- 


govofo 


2a 0 (l + ao) L d—1 


(d — 2 + Ao) 


P a p{p) + 


Mo ~ l ) 2 

d- 1 


-Tl a fl(3 


C d - i(p • n) 2 x 


-( 


£ 

(4.20) 


The renormalization constants are found from the re¬ 
quirement that the function (14.201) . when expressed in 
new renormalized variables, be UV finite, i.e., finite at 
£ —>• 0. From the analysis of this expression it follows, 
however, that the pole in £ in the structure with n a fip 


cannot be removed by renormalization of the model pa¬ 
rameters because the bare part of does not contain 
analogous term. In order to ensure multiplicative renor- 
malizability one has to add such term, with a new positive 
amplitude factor uq, to the bare part: 


-i a/3 


= {iui - z/ 0 pi - izq/o(p ' n) 2 } P a p(p) - v 0 f 0 u 0 (p • n ) 2 n a h 0 


go^ofo 


2a 0 (l + ao) |_ d—1 


(d — 2 + Aq) 


P a p( P) 


Mo ^ i ) 2 

d—l 


^ 0-^/3 


Cd- i(p • n) 2 x 


T~t 


(4.21) 


This means that the original model ED is extended by 
adding a new term of the form uofoi'o(n k 9 k )d 2 ,(n k 9 k )', 
the interpretation of the new parameter uq is literally 
the same as for / 0 in Sec. [TTJ 

Now the model is multiplicatively renormalizable with 
two independent renormalization constants Zf and Z u : 


u 0 = vZ v , / 0 = fZ f , u 0 = uZ u , 

A 0 = AZ a , g 0 = gn^ +ri Z g , a 0 = a/i 71 Z a , (4.22) 

I 


at that 

Z u = Z a = Z A = 1, Z 9 = Zj 1 * * * * . (4.23) 

Here g is the “reference mass” (additional free parameter 

of the renormalized theory) in the MS renormalization 
scheme, which we always use in what follows; g , u , a, 

v, A and / are renormalized analogs of the bare parame¬ 
ters go, Mo, ao, voi Ao and /o, and Zj = Zj(<7, £, d) are the 
renormalization constants. Their relations in (14.231) re¬ 
sult from the absence of renormalization of the contribu¬ 
tion with D~ x in (13.11) . so that D 0 = go^ofo = g// + V 3 /, 
a o^o = ag v u. No renormalization of the fields and the 
parameter mo = m is needed: i.e., Z$ = 1 for all and 

Z m = 1- 

The renormalized action functional has the form 


Sfl(4>) = \e[D e e' k - \ Vi d- 1 


Vk 


■0L 


-d t e 


t u k 


+ v fZf uZ u (n k 6' k )dl (n k 9 k ), 


{vidi)0 k + A(9idi)v k + v{d 2 1 _ + fZ f d^)6 k 


+ 


(4.24) 


where the function D v from (12.71) should be expressed in 
renormalized variables using (14.221) . 

At this moment one important problem springs up. 
Since the original model is extended by introducing a 
new term (proportional to the 9[9 k ) in the action func¬ 
tional ED, one may guess that the propagator func¬ 
tions ED and ED have to be modified. Conse¬ 
quently, we have to recalculate the diagrams for functions 
(9' a 9@) i-ir and (9' a 9pv 1 )i.i I , i.e., the expressions (14.131) 
and (14.171) . 

If fact, the difference between the original expressions 
for the bare propagators and the new ones is that the 
second have additional terms, which are proportional to 
the p||. Consequently, they do not contribute to the in¬ 


tegrals and revision of the final expressions is in fact not 
needed; this means that expressions (14.131) and (14.171) 
remain valid in the modified model. This problem was 
examined in details in [30]; moreover, the derivation of 
the propagators in the presence of a distinguished direc¬ 
tion n, i.e., in fact, the matrix inversion in the orthogonal 
subspace, was also discussed there. 

Now we are ready to study the fixed points {g*} that 
govern the IR asymptotic behavior. The basic RG equa¬ 
tion for a multiplicatively renormalizable quantity (cor¬ 
relation function, composite operator, etc.) has the form 

[Prg + 7-f] Fr = ® 


(4.25) 
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and is a consequence of operating on the relation F = 
Z f F r with the differential operation for fixed set of 
bare parameters eo = {go, Uh /o, uo, -To}- This operation 
is customarily denoted as V^, and 7 ^ is the anomalous 
dimension of F. Since Z v = \. the renormalization group 
operator T> R g has the form V R g = TV + P g d g — 7/2?/ — 
7 U TV, where V x = xd x for any variable x. 

The RG functions are defined as 


to 

III 

II 

to 

T 

1 

03 

1 

to 

S' 

(4.26a) 

Pu = VgU = -uiu(g,u), 

(4.26b) 

Pa = = -T)a, 

(4.26c) 


7 f = TV hr Z F = /3gd g In Zp for any Z F . (4.26d) 


The relations between p and 7 in (14.26al) (14.26cl) re¬ 

sult from their definitions along with relations (14.221) 
and (14.23[) . 

The constants Zi are found from the requirement of 
UV finiteness of the expression (14.211) . Thus, for the pa¬ 
rameter /o that splits the Laplace operator we obtain 


z f = \- 


{d - 2 + A) 
2(d — 1 ) 


g 1 
a(a + 1 ) £ 


+ o(g 2 ), 


(4.27) 


_ {d — 2 + A) g 
7/ “ 2{d — 1 ) a{a + 1 )’ 


(4.28) 


where we passed to the new coupling constant g = g Cd -1 
with Cd -1 from (14.131) . 

Then we have to renormalize the constant uq such that 
the expression 


go f quo 


1 + 


M-l ) 2 


1 




2(d — 1 ) uq ao(l + cko) £ 
be UV finite to the first order in g. Therefore, 
(A — l) 2 


n a np( p • n) 2 
(4.29) 


Z u Zj = 1 — 
and 


g 


1 


2 (d — 1 ) ua( 1 + a) £ 
(A-l) 2 


7 + 0(g 2 ), (4.30) 


7u + 7/ = 


2 (d — 1 ) u a{ 1 + a) ’ 


(4.31) 


where the constant 7 f is obtained in (14.281) . Furthermore, 
from the last relation in (14.231) it follows that for the 
coupling constant g 


(d ~ 2 + A) g 
2 (d — 1) a(l + a) 


(4.32) 


We stress that, since the expression (14.211) is exact, 
i.e., it has no corrections in coupling constant g , all the 
above expressions for the anomalous dimensions 7 f g u 
are exact, too. 


E. Fixed points 


One of the basic RG statements is that the asymptotic 
behavior of the model is governed by the fixed points 
{g*, a*, u*, /*}, defined by the relations 

Pg=0, fiu = 0 , 0f = 0 and fd a = 0 ; (4.33) 


here 


Pg = g (-£ 


V + 7/) = g 


-£-v + 


(d-2 + A) 
2(d — 1 ) 


g 

a(l + a) 
(4.34a) 


Pu = = 


9 

a(a + 1 ) 


(d-2 + A) V4-1) 2 1 

- —'ll — - 

2(d-l) 2(d — 1) J ’ 

(4.34b) 


Pf 


= _ M-2 + A) _ g _ 

Jlf J 2(d — 1 ) a(l + a) ’ 


(4.34c) 


the expression for /3 a is written in (I4.26cl) . 

The type of a fixed point (IR/UV attractive or a saddle 
point), i.e., the character of the RG flow in vicinity of the 
point, is determined by the matrix fVfc = dPi/dgk, where 
Pi is the full set of /3-functions and gk is the full set of 
couplings. For an IR attractive fixed point the matrix fl 
are positive, i.e., the real parts of all its eigenvalues are 
positive. 

The analysis of the /3-functions reveals several fixed 
points. The first possibility is to put a* = 0; conse¬ 
quently we get at once the trivial case g* = 0. There 
is, however, another possibility - to disclose it we have 
to pass from the coupling constant g to new constant 
g' = g/a , which is assumed to be finite at a —> 0. In 
fact this means, that the correlation function D v (co) be¬ 
comes proportional to <5(w) (see (12.71) 1 and we deal with 
the independent of time (“frozen” or “quenched”) veloc¬ 
ity field. 

The new /3-function, which remains nonzero at a — > 0, 
is 


Pg> - Jg 




(d — 2 + * 4 .) , 
2{d — 1 ) 9 


(4.35) 


the matrix fl in these variables has the form 


(9g'Pg> 

dg'Pu 

0 

dg'Pf 

0 

duPu 

0 

0 

0 

0 

da Pa 

0 

V 0 

0 

0 

dfPf 


This situation implies two options: 

(la) g'* = 0, with fl* g , g , = dp g >/dg'\ gl=g „ = -£ and 

Ka = -V 

For the two remaining parameters u and / we have 
Pu = Pf = 0, fluu = ^// = 0, so that both u and / 
remain free parameters. 

Since U*, u = 0, the matrix fl is triangle and its eigen¬ 
values coincide with the diagonal elements. Thus, this 
fixed point is IR attractive for £ < 0, 7 ? < 0; 
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(lb) if S'* = = £ and F!* Q = - 77 , so 

that this fixed point is IR attractive for £ > 0, 77 < 0. 
For the remaining parameters u and / we have the fixed- 
point values u* = {A — l) 2 / (d — 2 + A) and f* = 00 with 

Another interesting case to be considered is a* = 00 . 
From (12.7p it follows that this case corresponds to the 
rapid-change model with new charge g” = g/a 2 , which 
is supposed to be finite at a —> 00 . Besides that it is 
convenient to pass from the variable a to variable x = 
1/a, i.e., x —> 0. So, the new /3-functions are 

Px = xrj; (4.37a) 


= 9 


(d-2 + A) 
2(d — 1 ) U 


(A - 1) 2 ~ 

2{d — 1) J ’ 


(4.37b) 


Pf = 9" 


(d — 2 + Ad) 
* 2 (d-l) 


(4.37c) 



Pg" = ~o.P‘. 


a 


2 f 9 


29 ft " 
-jPa = 9 
a' 3 


-C + S + 


(d - 2 + A) „ 
2 (d-l) 9 ' 

(4.37d) 


Thus, we find two more fixed points: 

(2a) g"* = 0, with Q,*„ g „ = -£ + f?, Vt* xx = 77 . As in the 
case (la) for two remaining parameters u and / we have 
Pu = Pf = 0 , : f Yff = 0 , so both of them remain 

free parameters. 

As before the matrix Q in the new variables 
{g", x , u, /} is a matrix of the type (14.361) . i.e., it is trian¬ 
gle and its eigenvalues are simply given by diagonal ele¬ 
ments. Thus, this fixed point is IR attractive for 77 > 0 , 
il - £ > 0 ; 

( 2 b) if g"* = (f-77)^, n* g „ g „ = £-77 and Q* xx = 77 , 
so that this fixed point is IR attractive for 77 > 0, £ — 77 > 
0. For the remaining parameters u and / we have the 
fixed-point values u* = {A— l) 2 /(d — 2 + A) and f* = 00 

with n* uu = n* f = ^-ri. 

For the special case 77 = 0 the function p a and the 
eigenvalue Q aa vanish identically, so that the nontriv¬ 
ial fixed point [g/a(a + 1 )]* = 2 £(d — l)/(d — 2 + A) is 
IR attractive for £ > 0. Moreover, this fixed point is 
degenerate in the sense that we can not determine the 
parameters g* and a* separately. 

Thus, we can conclude, that the domains of IR stability 
in this vector model JS3D coincide with the correspond¬ 
ing domains of IR stability in scalar model, considered 
in [29|. The general pattern of the fixed points stability 
in the £ — 77 plane is shown in Fig. [T| The straight lines 
77 = 0 ; £ = 0 , 77 < 0 ; and £ = 77, 77 > 0 corresponds to the 
boundaries of domains, which has neither gaps nor over¬ 
laps between them. Since the /3-functions (14.341) have no 
higher-order corrections, this pattern is exact. 

Note that the Kolmogorov values of the exponents £ = 
8/3 ,77 = 4/3 lie deep inside the domain of stability of the 


FIG. 7. Domains of IR stability of the fixed points in the 
model m- The numbers in boxes correspond to the fixed 
points (la) - (26) in the text. 


nontrivial rapid-change point (26); there is no borderline 
going through this point. 

This fact implies that the correlation functions of the 
model (ED in the IR region (/rr ~ Ar !?> 1, Mr ~ 1) 
exhibit scaling behavior (as we will see below, up to log¬ 
arithmic factors). 

The corresponding critical dimensions A[.F] = for 
all basic fields and parameters can be calculated exactly; 
see the next subsection. 


F. Critical dimensions 

In the leading order of the IR asymptotic behavior the 
Green functions satisfy the RG equation (|4.25p with the 
substitution g —> g*, a —> a*, f —> f* and u u*. The 
operator is invariant with respect to the change of 
variables {x,y} ->■ {ai.e., p x d x + P y d y = P x ’d x f + 
Py’dy’. Taking into account the fact that 7 * = 0, this 
gives 


[V, - 7 }V f + 7 £] G R (e, /x, •••) = 0. (4.38) 

Canonical scale invariance is expressed by the relations 




G H = 0, 


Y. d ^- d C 


G H = 0, 


a 


a 


(4.39) 

where a = {t, x, g, v , a, m, M, u, /, A , g} is the set of all 
arguments of G R (t, x is the set of all times and coordi¬ 
nates), and d k and d u are the canonical dimensions of G R 
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and <7. Substitution of the needed dimensions from Ta¬ 
ble U and combination of the obtained result with (14.381) 
gives the desired equation of critical IR scaling for the 
model: 


and distinguishes it from both the isotropic Kraichnan’s 
vector model 0 (in which 7 ^ ^ 0 ) and anisotropic 
Kraichnan’s scalar model [29] (in which the Laplacian 
splitting parameter fo is not dimensionless). 


[—2? x + A ( P t + A m T> m + Am^m + A/l?/ — Aq] G r — 0, 

(4.40) 

where 

A t = -A*, = - 2 , A m = A m = 1 , 

A/ = 7/; A u = 0 (4.41) 

and 

A[G] = A G = 4 + 2 dg + 7 * G (4.42) 

are the corresponding critical dimensions. Substituting 
the values of fixed point of the regimes (la)-( 2 b) we ob¬ 
tain: 

A/=0 for (la), (2a); (4.43) 

A/ = £ for (lb), and A/ = £ — 77 for ( 2 b). 

In particular, for any correlation function G R = 
($...$) of the fields $ we have Aq = A$, with the 

summation over all fields < E > entering into G R , namely, 

Ag = _/V$d$ = Ng'de' + Ngdg + N v d v . (4.44) 

$ 

Since in the model m the fields themselves are not 
renormalized (i.e., 7 $ = 0 for all $, see sec. IlVDl) . us¬ 
ing (14.42[) we conclude, that the critical dimensions of 
the fields <I> = {v, 0 , 6'} are the same as their canonical 
dimensions, presented in the Table |TJ Namely, 

A„ = 1, Ag = -1, Ag' = d + 1. (4.45) 

It is the specific feature of the present model, which 
makes it similar to the zero-correlation time model [30] 


V. RENORMALIZATION AND CRITICAL 
DIMENSIONS OF COMPOSITE OPERATORS 

The analysis of the renormalization of composite oper¬ 
ators is nearly the same as in the rapid-change model [30|, 
so we will discuss it here very briefly. 


A. General scheme 

The central role in the following will be played by com¬ 
posite fields (“operators”) built solely of the basic fields 
6 : 

F Np = (6A) P ( n s 9 s ) 2m , (5.1) 

where N = 2 (p + m) is the total number of fields 9 , 
entering the operator. 

As was pointed out in 0 , the operator counterterms 
to a certain Fpr P involve only operators of the form (15.11) 
with the same value of N. Besides that, all the corre¬ 
sponding diagrams diverge logarithmically and one can 
calculate them with all external frequencies and momenta 
set equal to zero. 

Let us denote the closed set of operators, which can 
mix to each other in renormalization, as F = {Tat p }. 
The renormalization matrix Zp = {Zpjp^Np'} for this set, 
given by the relation 

Fn p = ^2 z Np,Np'F$ p/ , (5.2) 

p f 

is determined by the requirement that the 1 -irreducible 
correlation function 


( f n p ( x )°( x i) ■ ■ ■ 6( x iv))i_ ir = 


— z Np , Np ,(F]Vp'(x)9(x 1 ) . . . 9(x-n)) 1 _- II — ^Np, Np'^Np' (x;xi, • ■ ■, X N ) (5-3) 

p' p' 


be UV finite in renormalized theory, i.e., it has no poles 
in £ when expressed in renormalized variables (14.221) . 
This is equivalent to the UV finiteness of the sum 
J2 P ' ZnI^Np^Np'^O), in which 

T Np ’{x-9) = — dx 1 ... / dx N r Np '(x-,xi,.. .,x N ) 


x 9{x\)... 0 (xn) (5-4) 

is a functional of the field 9(x). 

The contribution of a specific diagram into the func¬ 
tional T Np i in (15.41) for any composite operator Fn p ' is 












represented in the form 


form 
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IVp' = V aP ... 1 $};-" 9 a 9 b ..., ( 5 . 5 ) 

where V a p... is the vertex factor, !“»"■ is the “internal 
block” of the diagram with free vector indices, and the 
product 9 a 9 b ■ ■ ■ corresponds to external “tails.” 

According to the general rules of the universal dia¬ 
grammatic technique (see, e.g., El), for any compos¬ 
ite operator F(x) built of the fields 9 1 the vertex V a R... 
in (1531) with k > 0 attached lines corresponds to the 
vertex factor 

V£ p (x-, xi,...,Xk ) = 6 k F Np (x)/S 9 (x 1 ) ...S 9 (x k ). ( 5 . 6 ) 

The arguments x\ ... x k of the quantity (15.61) are con¬ 
tracted with the arguments of the upper 9 ends of the 
lines (99') o attached to the vertex. 


B. Exact result for the diagrams 

Now let us turn to the calculation of the internal block 
/“)(■■■ of the diagrams. The one-loop diagram is repre¬ 
sented in Fig. ©• 



Y$ = T4 a i(k) V Z]b (-k)P ai (k)P Po (k)n x n z 

= -A 2 n x P xa (k)n z P z p(k) k a k b , (5.7) 


where the letters i. j, x and z denote internal indices of 
the diagram itself. Then we have to integrate over 
the frequency and momentum with the factors like m 
and m , namely 


dk 


1 


1 


rab _ 

a/3 J (27r) d —iuj + v\ k^_ + iuj + uk"^ -(- vfk | 

, 5-d-(^+?7) 

x 2^(fc||) Dq -i- -2 Y$. (5.8) 


Since the expression 




aovok A 

contains the factor <5(fc||), we 
can perform all the calculations with the original propa¬ 
gators m and (1371) : see the discussion in Sec. IIV Dl 
Using the relation (14.91) for averaging over the angles 
and setting rj = 0 [see the discussion after (14.101) ]. we 
arrive at the following result: 


rab 

l a(3 


A 2 f f dk± 1 fcq 
2a(l + a)^ J ( 2n ) d ~ 1 k\ a 13 


A 2 f 1 

2a(a + 1) (d — 1) 


Pab( n) n a np g x 


m ^ 

IT 


(5.9) 


Contributions of all multiloop diagrams are equal to 
zero, see Sec. nvn The multiloop diagrams of the “sand 
clock” type, represented by products of simpler diagrams, 
contain only higher-order poles in £ and, in the MS 
scheme, do not contribute to the anomalous dimensions. 
Therefore the one-loop approximation (15.91) gives the ex¬ 
act answer. 


FIG. 8. The one-loop contribution to the generating func¬ 
tional (15.41) . 


C. Renormalization matrix and anomalous 
dimensions 


Once all the external frequencies and momenta are set 
to zero, the index structure of this diagram takes on the 

I 


Combining expressions (15.51) . (15.61) and the exact an¬ 
swer (15.91) . for the functional T^p from (15.41) we obtain 


r^vp oc c/-i r/, [Fn p ] x Ti a Ti/3 x P ab (n) x 9 a 9 b — 

= 2m(2m — 1) x Fn p +i + (2 p + 8 pm — 2m(2m — 1)) x Fn p + 

+ (4p(p- 1) - 2p- 8pm) x Fjvp-i - 4p(p - 1) x F Np _ 2 , (5.10) 


up to an overall scalar factor. 


Expression (15.101) shows that the operators Fm p indeed 
mix in renormalization: the UV finite renormalized op- 
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erator F R has the form F R = F+ counterterms, where 
the contribution of the counterterms is a linear combina¬ 
tion of F itself and other unrenornralized operators with 
the same total number N of the fields, which are said to 
“admix” to F in renormalization. 

Let F = {.Fp} be a closed set of operators (15.11) with 
a certain fixed value of N (which we will omit below for 
brevity) and different values of p, which mix only to each 
other in renormalization. The renormalization matrix 
Zp = {Z p y} and the matrix of anomalous dimensions 
7 f = { 7 Pl p'} for this set are given by 

F p = J2 Z p yF R 7 7 p = Zp^Zp. (5.11) 

p' 

The scale invariance (14.391) and the RG equation (14.251) 
for the operator F p give the corresponding matrix of crit¬ 
ical dimensions A p = {A PjJ /} in the form similar to the 
expression (14.421) . where dp, dp and dp should be under¬ 
stood as the diagonal matrices of canonical dimensions 
of the operators in question (with the diagonal elements 
equal to sums of corresponding dimensions of all fields 
and derivatives constituting F) and 7 * = 7 (g*, a*,u*,f*) 
is the matrix (15.111) at the fixed point. 

In this notation and in the MS scheme the renormal¬ 
ization matrix Z has the form 


z = i + A, 


(5.12) 


where I is the unity matrix and the elements of the ma¬ 
trix A have the forms 


dpp/ — CL' 


VP 


(5.13) 


7 p. = 0. This means that for such £ and 77 the critical di¬ 
mensions of the composite operators coincide with their 
canonical dimensions, so that there is no corrections to 
ordinary scaling. 

For the regimes (lb) and (2b) we have g'* = 
and g"* = 73 ^ (£ - rj), so that 

7p,p'+i = 2 / X 2m(2m — 1); 

7 p, P ' =V x [ 2 p + 8 pm - 2m(2m - 1 )]; 

7 p, p'-i = V X [4 p{p - 1 ) - 2 p- 8 pm] ; 


7p,p '- 2 = V x [~Mp ~ !)], (5-16) 


where y denotes the common factor, i.e., 


y = - 


A 2 f 

2(d — 2 + 71) 




for the critical regime (lb); 

(5.17a) 


A 2 f 

V = 2 (d-^'+A) ^ ~ ^ f ° r the critical re § ime ( 2b )- 

(5.17b) 

Therefore the matrix of critical dimensions for the set 
F p with fixed N has the form 


Ap lP ' = -2(p + m)5 pp f + 7 * iP /, (5.18) 

where —2 (p + m) is the canonical dimension, 5 PP ' is Kro- 
necker’s 6 symbol and 7 * , is the value of the matrix of 
anomalous dimensions at the fixed point. 


D. Asymptotic behavior of the correlation function 

G = (F 1 F 2 ) 


Since the renormalization matrix Z has the 
form (15.121) , the matrix of anomalous dimensions 7 
has the form 


7 pp 1 — ctpp' g 


(5.14) 


with the coefficients a pp ' from (15.131) . Combining (15.101) 
(15.141) and taking into account the scalar factor, not writ¬ 
ten in (15.101) . but presented in (15.91) . together with the 
fact, that the symmetrical coefficient for this one-loop 
diagram is 1 / 2 , one obtains the following expression for 
the matrix of anomalous dimensions 7 : 


7 p, p '+1 


7 p,p' 


7 p, p'-i 


7 p, p '—2 


A 2 f 


1 


Aa{a + 

1 ) 

(d- 

1 ) 

A 2 f 


1 


Aa(a + 

1 ) 

(d~ 

1 ) 

A 2 f 


1 


4 a(a + 

1 ) 

(d- 

1 ) 

A 2 f 


1 


Aa{a + 

1 ) 

(d- 

1 ) 


2m(2m — 1 ) g\ 

[ 2 p + 8 pm — 2m(2m — 1 )] g\ 
[4 pip - 1) - 2p - 8 pm] g ; 
[~4p(p-l)\ g. (5.15) 


Now we have to substitute the value of the fixed point 
into the expressions (15.151) . For the critical regimes 
(la) and ( 2 a) we immediately arrive at the trivial result 


Up to a scalar factor y , the values of the matrix el¬ 
ements of the matrix of anomalous dimensions at the 
fixed point (| 5.16 1) are the same as in the zero-time cor¬ 
relation case [30| |. This means that the matrix of critical 
dimensions (15.181) is not diagonalizable, but can only be 
brought to the Jordan form, i.e., A p = UpApUp 1 , where 
the matrix Ap is 


A p = 


(—2 (p + m) 1 0 

0 — 2(p + m) 1 

: 0 


V 


0 


0 —2(p + m)J 
(5.19) 

For the equal-time pair correlation function of two 
composite operators Fpr p of the form (15.11) with arbitrary 
values of N and p 


G n 

1 P 1,N 2 P2 (r) = (Fpr lPl (t,x 1 ) F n 

2P2 (A x 2 )), (5.20) 

where r = |x 2 — xi|, i = {IVipi}, and k = {N 2 P 2 }, this 
leads to the appearance of logarithmic dependence in the 
IR asymptotic behavior (in the following we denote in 
G ik for brevity): 
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Gf k oc (^) JVl+Ar2 -P(Ar 1 +Ar 2 )/2 [Innr] $ (l, Mr, mr, f) Vi, k. 


(5.21) 


Expression (15.211) is written up to a dimensional constant 
factor; Pl (...) is a polynomial of degree L with the ar¬ 
gument In nr', f is the invariant charge and / —>- fr^ 
as 1 /nr —t 0 for scaling regime (lb), / —>• fr^~ v as 
l/n r ~> 0 for scaling regime (2b). 

Representations (15.211) with yet unknown scaling func¬ 
tions *3> (Mr, mr, /) = $ (l, Mr, mr, f) describe the be¬ 
havior of the correlation functions for /rr 1 and any 
fixed value of Mr. The inertial range £ <C r<Cl corre¬ 
sponds to the additional condition Mr <C 1. Here and 
below we do not distinguish the two IR scales M and m, 
first introduced in (12.21) and m-, the form of the func¬ 


tions $ (Mr, f) 


as Mr —>• 0 is studied using the 

/=const 


J 


operator product expansion. 

In general, the operators entering into the OPE are 
those which appear in the corresponding Taylor expan¬ 
sions, and also all possible operators that admix to them 
in renormalization In our case the main con¬ 

tribution to the sum is given by the operator F R oc 
(Mr) Nl+N2 x P(n 1 +n 2 )/2 [In Mr] which possesses max¬ 
imal singularity. 

Combining this fact with the RG representation (15.211) , 
restoring canonical dimension da = — N\ — N 2 and re¬ 
taining only the leading term, we obtain the following 
asymptotic expression for the pair correlation function 
G 1)5.201) in the inertial range: 


G = (Fn, Fn 2 


OC /gM-'Vi-V [l n nr ](N 1 +N 2 )/2 \\ nMr \(Ni+N 2 )/2 


(5.22) 


where (/) is a certain scaling function, restricted in 
the inertial range f « r < L. Owing to the nilpotency 
of the matrix of critical dimensions, the result obtained 
is independent of the scalar factor y (15.171) , and the only 
dependence on the exponents £ and rj, that distinguishes 
two nontrivial cases (lb) and (2b), is contained in the 
invariant charge /. 

For the trivial regimes (la) and (2a) there is no cor¬ 
rections to ordinary scaling. 

VI. CONCLUSION 

We applied the field theoretic renormalization group 
and the operator product expansion to the analysis of the 
inertial-range asymptotic behavior of a divergence-free 
vector field, passively advected by strongly anisotropic 
turbulent flow. 

Depending on the two exponents £ and rj that describe 
the energy spectrum £ oc kp ^ and the dispersion law 
w ~ kp 71 of the velocity field, the possible nontrivial 
types of the IR behavior appear to reduce to only two 
limiting cases: the rapid-change type behavior, realized 
for £ > rj > 0, and the “frozen” (time-independent or 
“quenched”) behavior, realized for £ > 0, rj < 0. 

To avoid possible confusion we stress that we studied 
the model with arbitrary finite correlation time of the 
velocity field. The behavior typical of the vanishing or 
infinite correlation time is formed effectively in the IR 
range as the leading-order asymptotic behavior of the 
correlation functions. 

In this respect, the situation is the same as in the 
model of the anisotropic advection of the scalar field, 
studied in [29). Thus, another important conclusion 


of that work remains true - in contrast to the frnite- 
correlated isotropic case, where the Kolmogorov values 
£/2 = y = 4/3 lie exactly on the crossover line between 
the rapid-change and frozen regimes [Hj, 0 E} , in the 
present model they lie inside the domain of the rapid- 
change regime; there is no crossover line going through 
this point. This result is in agreement with the analy¬ 
sis of [13 and in disagreement with the [U, [25| for the 
scalar case. 

The inertial-range asymptotic expressions for vari¬ 
ous correlation functions are summarized in expres¬ 
sions (15.221) . In contrast to the Kraichnan’s rapid-change 
model, where the correlation functions exhibit anoma¬ 
lous scaling behavior with infinite sets of anomalous 
exponents, here the dependence on the integral turbu¬ 
lence scale L demonstrates a logarithmic character: the 
anomalies manifest themselves as polynomials of loga¬ 
rithms of ( L/r ), where r is the separation. 

The key point is that the matrices of scaling dimensions 
of the relevant families of composite fields (operators) 
appear nilpotent and cannot be diagonalized - they can 
only be brought to Jordan form; hence the logarithms. 
This result is perturbatively exact in the sense that the 
contributions of all multiloop diagrams appear equal to 
zero. 

The possibility of logarithmic dependence of various 
correlation functions on the integral scale L and the sepa¬ 
ration r should be taken into account in analysis of exper¬ 
imental data. Since the difference between the nontrivial 
regimes (lb) and (2b) stays only in the argument of the 
scaling function $, it requires very accurate experiments 
to discern them. 

It remains to admit that, although our model has a 
finite correlation time and possess Galilean symmetry, it 
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is still simplified in the sense that the velocity ensemble 
is Gaussian. More realistic models should involve the 
nonlinear NS equation, while the anisotropy should be 
introduced by the large-scale stirring. So far, the analysis 
based on the advecting NS velocity field was performed 
only for the passive scalar [Hj and vector [Hj fields only 
in isotropic cases. 

Thus, the analysis of the full-scale problem remains for 
the future; this work is already in progress. 
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